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ABSTRACT 

We reanalyzed the ROSAT/HRI observation of MS1054-03, optimizing the channel HRI selection and 
including a new exposure of 68 ksec. From a wavelet analysis of the HRI image we identify the main 
cluster component and find evidence for substructure in the west, which might either be a group of 
galaxies falling onto the cluster or a foreground source. 

Our 1-D and 2-D analysis of the data show that the cluster can be fitted well by a classical /3-model 
centered only 20" away from the central cD galaxy. The core radius and /3 values derived from the 
spherical model(/3 = 0.96±g;||) and the elliptical model (/3 = 0.73 ± 0.18) are consistent. 

We derived the gas mass and total mass of the cluster from the /3-model fit and the previously published 
ASCA temperature (12. keV). The gas mass fraction at the virial radius is f gas = (14[— 3, +2.5] ± 
3 )% for Qq — 1, where the errors in brackets come from the uncertainty on the temperature and the 
remaining errors from the HRI imaging data. The gas mass fraction computed for the best fit ASCA 
temperature is significantly lower than found for nearby hot clusters, f gas = 20.1 ± 1.6 %. This local 
value can be matched if the actual virial temperature of MS1054-032 were close to the lower ASCA limit 
(~ 10 keV) with an even lower value of 8 keV giving the best agreement. Such a bias between the 
virial and measured temperature could be due to the presence of shock waves in the intracluster medium 
stemming from recent mergers. Another possibility, that reconciles a high temperature with the local 
gas mass fraction, is the existence of a non zero cosmological constant. 

Subject headings: galaxies: clusters: general, individual (MS1054-03) — X-rays: galaxies — cosmology: 
dark matter and observations 



1. INTRODUCTION 

MS 1054-03 is the most distant cluster of galaxies found 
in the X-ray selected sample of the Einstein Extended 
Medium Sensitivity Survey (EMSS; Gioia et al., 1990). 
Due to the extreme sensitivity of the high-mass end of 
the cluster mass function to the density parameter, Hq, 
the existence of massive, virialized clusters at high z is of 
great cosmological significance (e.g. Oukbir & Blanchard 
1992). The very existence of even a few massive clusters 
at redshifts approaching unity and beyond strongly argue 
against cosmological models with f2o = 1. 

There is a general consensus that MS1054-03 is such 
a massive cluster, based on its high X-ray temperature, 
kT = 12.3lJ:2 keV > measured with ASCA (Donahue et 
al. 1998 hereafter D98), its high velocity dispersion (Tran 
et al. 1999) and intense weak lensing signal (Luppino & 
Kaiser 1997). On the other hand its apparent complex 
morphology (D98) might indicate that MS1054-03 has not 
yet reached an equilibrium state, which could cast some 
doubt on the dynamical mass estimates and on the rel- 
evance of the cluster for cosmological tests. Tran et al. 
(1999) emphasized the good agreement between the vari- 
ous mass estimates, but the large error bars must be noted. 

In this paper we re-investigate the ROSAT/HRI ob- 
servations of MS1054, including a new exposure taken in 
1997. We first try to better understand the cluster mor- 
phology. A wavelet analysis of the HRI image is performed 
to unravel significant substructures and identify the main 
cluster component. A /3-model is then fitted to the data. 
From the fit results, we estimate for the first time the gas 



mass of this cluster. A comparison of the gas mass frac- 
tion of MS 1054-032, with the gas mass fraction of nearby 
clusters, provides a consistency check on the total mass 
estimate, assuming that the gas mass fraction in clusters 
does not evolve with redshift. 

Throughout the paper we assume Hq = 50 km/sec/Mpc, 
qo = 0.5, A = and all quoted error bars are 1 a unless 
otherwise stated. At the redshift of the cluster, z — 0.83, 
one arcmin corresponds to 497h 50 1 kpc. 

2. OBSERVATIONS 

MS1054-03 was observed by the ROSAT HRI (Triimper 
1992) for 190,754 sec in total: an exposure of ~ 122 ksec 
in 1996, which was analyzed and presented in D98 and 
a pointing of 68 ksec taken in 1997. In this analysis we 
combined the exposures of 1996 and 1997 and we selected 
only HRI channels 1 to 7 (David et al. 1997) in order to 
maximize the signal-to-noise ratio. The background level 
is thus lowered by a factor of about 2. The HRI image 
smoothed with a Gauss filter is shown on Fig. 1 and looks 
similar to the image presented by D98. 

3. MORPHOLOGY 

In order to remove noise and to identify the significant 
structures we performed a wavelet-analysis of the HRI raw 
image with the Multi-scale Vision Model (MVM) package 
(Rue & Bijaoui 1997). We explicitly took into account 
the Poisson statistics and the significance level was set to 
3o\ The algorithm is optimized for the detection of dif- 
fuse components by normalizing the wavelet coefficients 
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by their energy (Rue & Bijaoui 1997, Arnaud et al. 2000). 

The reconstructed image consists of point sources and 
one large scale diffuse component at the cluster posi- 
tion. To extract possible substructures, we reapplied the 
multi-scale analysis to the reconstructed diffuse source 
only. Two components are detected that are shown in 
Fig. 2, superimposed on the cluster optical image. The 
main component (~ 90% of the flux) is centered at RA 
= 10 /l 56 m 58.47 s , Dec.= -03deg37'31.2", about 20 " 
(166 h^, 1 kpc) from the dominant galaxy (D98). We iden- 
tify this structure with the main cluster component. It 
is very elliptical with an elongation in the east-west di- 
rection. The orientation follows the filamentary structure 
found in the optical and in the weak lensing map (Luppino 
& Kaiser 1997). The second component is compact and 
lies in the west, with an offset of 0.6 arcmin with respect 
to the main cluster center. It coincides with the brightest 
image peak, identified by D98, but its reconstructed flux 
is only ~ 10% of the main cluster flux. 

This western source might be a subcluster gravitation- 
ally connected to the main cluster, or a projection effect 
due to a group of galaxies in the field-of-view (FOV) or 
even a point source. The available data are insufficient to 
settle this issue. There is a weak indication that the source 
is extended, from our comparison of the source surface 
brightness profile with the point spread function (PSF) of 
the ROSAT/HRI (David et al. 1997). However, due to 
the lack of bright sources in the FOV, we cannot correct 
for possible alignment uncertainties, and the extent might 
be due to errors in the aspect reconstruction or residual 
main cluster emission. Very close to the maximum of this 
source we find several cluster members so this structure 
might be a subgroup falling right onto the cluster. Finally 
it must be noticed that the reconstructed surface mass den- 
sity by Luppino & Kaiser (1997) presents an extension in 
the west. Although it is unresolved, it might indicate that 
the western X-ray component does trace a gravitational 
potential. 

4. ISOTHERMAL BETA-MODEL FITS 

We performed /3-model fits to the data using spherically 
symmetric (1-D) and elliptical (2-D) /3-model fits. We 
excluded from the fits circular regions corresponding to 
serendipitous point sources. In the 1-D case the X-ray 
surface brightness profile is given by: 

S = S (l + r 2 /rty 3fj+ °- 5 + B (1) 

where r c is the core radius, /3 is a slope parameter and 
B is the background level. The functional form for the 
elliptical /3-model can be found in Neumann & Bohringer 
(1997). 

4.1. Spherical model 

We binned the data in concentric annuli with the center 
at the position determined from the above wavelet recon- 
struction. The /3-model was convolved with the PSF of the 
ROSAT/HRI. The free parameters in the fit are So, r c , (3 
and B. We first consider the main cluster component in a 
common way by excluding non-symmetric features in the 
data, i.e. a circular region around the substructure in the 
west. Fig. 3 shows the observed surface brightness profile 
and the best fit model. The reduced % 2 value is \ 2 = 0.97, 



and the shape parameters are well constrained (see Tab .[!]). 
On the other hand if the substructure is not removed, a 
good fit is still obtained (x 2 — 0.99) but the best fit shape 
parameters are unreasonably large ((3 = 2.,r c = 840 kpc), 
and poorly constrained (/3 > 1.1). An artificial increase 
of the best fit shape parameters is a well known effect of 
not excising a secondary sub-cluster in the /3-model fit 
(Jones & Forman 1999). A good fit is still obtained due to 
the large errors on the observed profile, whereas no upper 
limit on /3 can be set because the best fit core radius be- 
comes similar to the maximum radius of cluster emission 
detected. 

4.2. Elliptical model 

We also fitted an elliptical isothermal /3-model to the 
HRI image. Again we excise the substructure in the west 
as well as other weak unresovled sources in the FOV. Due 
to the limited statistics of the observation, we did not at- 
tempt to fit the whole image with a /3-model and an addi- 
tional component for the substructure, since the elliptical 
/3-model has already 8 fit parameters. 

The iso-contours of the best fit model are superimposed 
on the cluster image in Fig. 1. The fitting procedure is 
described in full detail in Neumann (1999). We binned 
the data into an image with a pixel size of 5" x 5". The 
fit included all pixels less than 4.6 arcmin away from 
the central pointing position of 10' l 57 m 00.00 s , and Dec.= 
-03deg37T2.0" (J2000). The central position of the clus- 
ter was left free to vary. 

As we are dealing with low number statistics in each 
image pixel, which shows non-Gaussian behavior, we 
smoothed the data with a Gauss filter (a = 10") before 
fitting. The modeling accounts for the Gauss filtering and 
for the PSF (the fitted model is convolved with the Gauss 
filtered PSF). To calculate the errors of the /3-model pa- 
rameters we performed a Monte-Carlo analysis in which 
we added Poisson noise to the data and fitted the /3-model 
subsequently. We performed 100 Monte-Carlo realizations. 

The validity of this approach was discussed in Neumann 
& Bohringer (1997) and Neumann (1999). In order to see 
whether it is still free of systematics in this regime of ex- 
tremely low signal-to-noise data (the central cluster inten- 
sity is only a factor of two higher than the background 
level) we simulated 100 realizations of the image of the 
best fit /3-model cluster including background. The Pois- 
son statistics were defined according to the length of the 
actual exposure time. We subsequently fitted a /3-model 
to these artificial smoothed images, as for the real image. 
The results are satisfactory, as the mean output values 
are practically identical to the input values, with differ- 
ences much smaller than the actual determined error bars 
for MS1054-03. The width of their distribution is slightly 
smaller than the errors determined from the real data. 

4.3. Comparison of the 1-D and 2-D models 

The values of the spherical and elliptical /3-model pa- 
rameters for the main cluster component together with 
their 1 a errors are given in Tab.[l] and the two /3-models 
profiles are plotted on Fig. 3. The best fit center of the 
elliptical model is only 5 " away from the central position 
from the wavelet analysis, which we also choose as the cen- 
ter for the 1-D fit. As already mentioned, these central 
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positions are also in excellent agreement with the posi- 
tion of the brightest cluster member. The 1-D and 2-D 
best fit parameters are consistent. The background level 
is slightly higher in the 1-D fit which explains the lower 
central value and the higher value for (3. Nevertheless the 
68% confidence intervals for the fitted shape parameters 
show a large large region of overlap. 

5. LUMINOSITY 

The background estimated from the 1-D model fit was 
subtracted from the radial profile. An additional 5% sys- 
tematic uncertainty was assumed on its level. The cluster 
emission is detected up to 2' or 1 h^, 1 Mpc at the 68% con- 
fidence level. The observed count rate within this aperture 
is (9±2) x 10 -3 cts/s, ~ 25% higher than the value given by 
D98 (subtracting additional sources, including the western 
structure for consistency). We recall that the count rate 
of the western substructure is only 10% of the cluster flux. 
The count rates derived from the best fit 2-D model is 
consistent but higher 11 x 10 _3 cts/s, a consequence of 
the lower best fit background, in this case. 

The observed count rate was translated into luminos- 
ity, assuming a cluster temperature of 12.3 keV, as mea- 
sured by D98 from ASCA data. The derived X-ray lu- 
minosity is L x = 1.2 x 10 45 h5 2 crgs/s (0.1-2.4 keV rest 
frame). The bolometric luminosity within 2' in radius, 
Lho\ = 4.3 ± 0.9 x 10 45 h^" 2 ergs/s, is consistent with the 
ASCA value of 4.4 x 10 45 h^ 2 ergs/s (D98). This might 
be fortuitous though since the later includes all compo- 
nents in the field. The bolometric luminosity is lower than 
the value expected from the Lx~T relation of Arnaud & 
Evrard (1999), assuming this relation does not evolve with 
rcdshift. For a cluster with T = 12. 3 jl§' 2 keV this correla- 
tion predicts a luminosity of 9.1T.|;o x 10 45 h^~ 2 ergs/s. Up 
to now it is not clear whether the Ly^-T relation evolves 
with rcdshift or not (see for example Schindlcr 1999). Also, 
the total X-ray luminosity of MS1054 might be actually 
higher since the cluster might extend beyond the radius of 
detection. However, due to the low S/N ratio of this obser- 
vation such extrapolation is very sensitive to the assumed 
background level and /3 parameter. 

6. MASS CONTENT 

6.1. Mass estimates 

The cluster gas mass profile M gas (r) can be derived from 
the best fit /3-model , given the observed temperature and 
Nh values. As the emissivity in the ROSAT/HRI energy 
band is nearly insensitive to the temperature, the uncer- 
tainty on the gas mass is overwhelmingly dominated by 
the errors on the /3-model parameters. The association of 
the western source with the cluster is unclear. However, 
in practice, this additional uncertainty has no significant 
impact on the gas mass estimate. If the substructure is 
included in the /3-model analysis, the derived gas mass 
differs by less than 5% from the value obtained when ex- 
cising it; the derived mass is increased within the detec- 
tion radius (as a consequence of the ~ 10% higher flux) 
and artificially decreased when one extrapolates the data 
to higher radii (due to the higher derived (3 value). In the 
following we only consider the main cluster component and 
its corresponding /3 model. 



The total mass can be estimated with the isothermal 
/3-modcl approach (BM), which is thought to be roughly 
valid even if the cluster is not fully in hydrostatic equilib- 
rium (e.g Schindler 1996). The alternative is to employ 
the virial theorem (VT) at given density contrast, over 
the mean mass density of the Universe at the cluster red- 
shift, normalized from numerical simulations (Evrard et 
al. 1996). The resulting M-T relation (see Appendix) de- 
pends both on redshift and on the cosmological parameters 
Slo and A (Voit & Donahue 1998). We used the analytical 
expression from Bryan & Norman (1998). 

We first fix the temperature to the best fit ASCA value 
of kT = 12.3 keV and consider an Qq = 1 Universe. The 
uncertainty introduced by the errors on the temperature 
is discussed later. The statistical errors on M gas (r) and 
on MbmW due to the uncertainties in the /3-model , are 
estimated following the general method described in El- 
baz et al. (1995). The derived mass profiles are plotted in 
Fig. 4. Tab. 2 summarizes the mass estimates at 1 h^~ Mpc 
(the maximum radius of detection) and 1.65 h^ 1 Mpc (the 
virial radius Ry for z=0.83 and kT = 12.3 keV using the 
Evrard et al. (1996) simulations). 

The 2-D model gas mass at 1 h^, 1 Mpc is 15% higher 
than the value derived from the 1-D model. This dis- 
crepancy, larger than the formal statistical uncertainties, 
is essentially due to the systematic difference in the back- 
ground estimates, which dominates the error. As the de- 
rived /3 value is consistently smaller in the 2-D fit, this 
discrepancy is amplified for extrapolated gas masses. It 
reaches 25% at the virial radius. Similarly the total BM 
mass estimate which scales as /3 is smaller for the 2-D 
model than for the 1-D model. Both are consistent with 
the VT estimate, which we will adopt in the following dis- 
cussion. The gas mass is taken as the average of the 1-D 
and 2-D estimates. Their difference is an estimate of the 
systematic uncertainties that we add in quadrature to sta- 
tistical ones. The gas mass fraction at the virial radius is 
thus f gas = 14 ± 3%, for kT = 12.3 keV and Q a = 1, the 
error coming from the uncertainty on the gas mass from 
the imaging data. 

We performed the same analysis for the extreme val- 
ues of the temperature, as allowed by the ASCA data 
at the 90% confidence level (D98). For the lower limit 
on the temperature, kT = 10.1 keV, the virial radius 
is 1.58 Mpc and the gas mass fraction within that ra- 
dius is f gas = 16.5 ± 3%. For kT = 15.4 keV, we get 
i?v = 1-95 Mpc and f gas = 11 ± 3%. 

In summary the gas mass fraction at the virial radius is 
fgas = 14[-3, +2.5] ± 3 % for fi = 1- We have separated 
the uncertainties due to the errors on the temperature (in 
bracket), which affect essentially the total mass estimate, 
and the uncertainties from the imaging data, which only 
affect the gas mass estimate. 

6.2. Comparison with nearby clusters 

As mentioned above, a precise determination of the total 
mass of distant luminous clusters like MS1054-03 is crucial 
for the determination of f2o based on cluster abundances 
at high redshifts. Unfortunately there are still important 
uncertainties on the total mass estimate: statistical errors 
due to the errors on the temperature measurement but 
also possible systematic errors if the temperature is a bi- 
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ased estimator of the cluster mass. In particular the gas 
temperature might differ from the true virial temperature 
if the cluster were not in hydrostatic equilibrium. 

By considering the additional information on the gas 
mass, we can however perform a consistency check on the 
total mass estimate. From a study of the intrinsic disper- 
sion in the Ly^-T relation, Arnaud & Evrard (1999) showed 
that the fractional variations of f gas at fixed cluster mass 
is very small. Since it is unlikely that the gas mass frac- 
tion evolves with redshift, the derived gas mass fraction of 
MS1054-03 should be consistent with the typical value of 
hot nearby clusters. 

As a reference we first consider the best fit ASCA tem- 
perature (kT = 12.3 kcV) and a fl = 1 universe. The 
corresponding gas mass fraction of MS1054-03 is signifi- 
cantly smaller (at the 95% confidence level) than the value 
found by Arnaud & Evrard (1999) for hot nearby clusters, 
f gas = 20.1±1.6 %. Obviously MS1054-03 might be a true 
outlier. Such outliers are rare but some, such as A1060, a 
cluster of exceptional low gas content in the local Universe 
(Arnaud & Evrard 1999), do exist. If this is not the case 
for MS 1054, the low derived f gas value suggests that the 
actual virial temperature is lower than 12.3 keV and/or 
that the adopted cosmology is wrong. We examine each 
possibility in turn. 

The derived gas mass fraction depends on the assumed 
values of the cosmological parameters fio and A. In the 

3/2 

BM approach, f gas oc g? a , where g?a is the angular dis- 
tance (Pen 1997). In the VT approach adopted here the 
dependence of f gas is slightly different. f gas is both de- 
pendent on the geometry factor (the variation of cZa) and 
on the normalization of the My(orRy)-T relation. A de- 
tailed derivation of this dependence and necessary equa- 
tions to compute the variation of f gas with Q are given 
in the Appendix. The variation of the gas mass frac- 
tion of MS1054-03 with tt is plotted on Fig. 5 for an 
open Universe ( Q < 1, A = 0) and a flat Universe 
(n + A = 1). For the best 1-D(2-D) /3-modcl , the esti- 
mated gas mass fraction of MS1054-03 would be a factor 
of 1.2(1.3) larger if fl = 0.3, A = 0.0 and 1.4(1.5) times 
larger if £lo = 0.3, A = 0.7. The gas mass fraction of 
MS1054-03 would in the later case be perfectly consistent 
with the local value, which would itself remain essentially 
unchanged (< 5% increase). 

The second possibility is that the virial temperature, 
and thus the virial mass, are smaller than given by the best 
fit ASCA value. The lower limit on the ASCA temperature 
(kT = 10.1 keV) yields a gas mass fraction marginally 
consistent with the local value. If the virial temperature 
were actually even somewhat lower, kT = 8 keV, the 
virial mass would decrease to 0.84 x 10 15 Mq, the virial 
radius to 1.3 Mpc and the gas mass within that radius, 
1.7 x 10 15 Mq, would reach 20% of the virial mass and 
be a perfect match with the local value. A virial tem- 
perature of 8 keV is formally excluded by ASCA mea- 
surements at the 90% confidence level but would be in 
better agreement with the measured velocity dispersion of 
cr = 1170 ± 150 kms/s (Tran ct al. 1999). In that case 
MS1054-03 would fit perfectly in the a — T relation estab- 
lished for z = 0.19 — 0.55 clusters and the virial relation 
(Fig. 3 of Tran et al. 1999). The temperature would also 



be more consistent with the measured bolometric lumi- 
nosity. A possible explanation for the measured tempera- 
ture being significantly higher than the virial temperature 
(apart from contamination by a hard source like an ab- 
sorbed AGN) would be the presence of shock waves in the 
gas if the cluster were indeed undergoing a recent merger. 

7. CONCLUSION 

Our wavelet analysis finds evidence for two components 
in the X-ray image of MS1054-03: a main diffuse com- 
ponent, with emission peaking at 20" from the brightest 
cluster member, and a compact substructure in the west, 
which coincides with the X-ray maximum in the ROSAT 
image. We emphasize that this brightest peak is not the 
centroid of the cluster. Indeed, when the western struc- 
ture is excluded, the cluster emission is fitted well by a 
classical /3-model , centered within 20" of the central cD 
galaxy. This indicates that we are then identifying the 
main cluster component. 

Unlike previous attempts (D98, Ebeling et al. 1999), we 
obtained a good fit to the data with a /3-model (low \ 2 
value) and derived well constrained parameters. This is a 
natural consequence of the higher S/N ratio of our data 
(optimum channel selection and longer exposure time) and 
our identification and removal of the western structure 
from our fits. 

There are indications from the X-ray data that the clus- 
ter is not fully relaxed. The core radius (~ 400 kpc) and 
cllipticity are relatively high and the western substruc- 
ture could be an in-falling group. This is not unusual, 
such mergers do exist in the nearby Universe. Actually 
MS1054-03 appears very similar to A521 at z=0.27 (Ar- 
naud et al. 2000), where the brightest peak is also associ- 
ated with the subcluster. 

The gas mass fraction derived for the best fit ASCA 
temperature of kT = 12.3 keV is only consistent with 
the local value if Qq < 1, a flat A-dominated Universe 
being favored. The local value can be matched as well 
for Clo = 1, provided that the actual virial temperature is 
close to the lower ASCA limit (~ 10 keV), with an even 
lower value of 8 keV giving the best match. To decide be- 
tween the two options requires a systematic analysis of a 
large sample of distant clusters. MS1054-03 would appear 
as an outlier in the second case (low virial temperature) 
and not in the first case (flat A dominated Universe). If 
the cluster's actual mass is indeed lower than previously 
estimated, this might have consequences for the measure- 
ments of f2 based on the abundance of massive clusters 
at high redshift. 

Finally we want to stress that the current data allow us 
to constrain the physics of MS1054 only with large error 
bars. Therefore better data, in particular spectroscopic 
data from Chandra and XMM, are strongly needed. 
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In this Appendix we consider the gas mass fraction derived from X-ray data in the VT approach (see Sect. 3.1) and 
examine how it depends on the assumed cosmological parameters, f2o and A. 

The total mass My is estimated from the measured temperature and the mass-temperature relation derived from the 
virial theorem at given density contrast: 

M v ocfM^»y 1/2 (1 + 2) -3/ 2(kT) a/ 2 (1) 



Q(z) 

It depends on f2 and A, via the density contrast, A c (f2(z), A), and the density parameter of the universe at redshift z, 
f2(z). The analytical expression of A c and A(z) can be found in Bryan & Norman (1998). The corresponding virial radius 
scales as: 



We assume that the gas density n g (r), follows a /3-model: n s (r) = n g ,o [l + (r/r c ) 2 ] 3 ^ 2 ■ X-ray imaging data provides 
the central surface brightness So, the slope parameter /3 and the angular core radius 8 C = r c /dx, where cLa is the angular 
distance. 

The emission measure along the line of sight through the cluster center can be derived from So, independently of any 
cosmological parameters : 

_ 4^(l + z) 4 S 

EM ° ~ A(T,z) (3) 

where A(T, z) is the emissivity in the detector band, taking into account the interstellar absorption and the instrument 
spectral response. 

Assuming that the X-ray atmosphere extends up to the virial radius, the central emission measure along the line of 
sight is linked to the gas density by EMq oc J q Rv n g (r) dr, whereas the gas mass within Ry is M gas oc J Q Rv n g (r) r 2 dr. 
By combining the two expressions, we can derive an expression for the gas mass that varies as: 

M gas ai? v /2 ^EM~ Q Q(fl, x c ) (4) 



where we have introduced the form factor Q(f3,x c ) = (n s )/y (w|); os . Here (n s ) is the average gas density within Ry 

and (n g )i os is the average along the line-of-sight passing through the cluster center. For the /3— model, this form factor 
depends on (3 and the scaled core radius x c = r c jRy\ 

n(n , *ti[l + {xlx c )*\- 3P ' 2 x 2 dx BETACFd.f,^) 

Q{P, x c ) = — " 1 == ^- = == (5) 



^^[l + ix/x^]-^ dx v /BETACF(i,3/3-i, T ^; 



BETACF is the continued fraction entering the expression of the incomplete Beta function B x (a,b) (see Press et al. 
1986) with: 

B x (a, b) = ~ X ^ BETACF(x, a, b) (6) 

From EqJ^, Eq.|, Eq.| and Eq.|, the gas mass fraction, estimated from given X-ray data, depends on the assumed 
cosmological parameters Slo and A as: 

fz &s aR\j 2 Q{p,e c d A /R Y ) (7) 

— 1 /2 

where Ry oc (A c (Q(z), A)£lo/Q(z)) 1 and g?a is the angular distance. Note that the variation with the cosmological 
parameters is not the same for all clusters. It depends on the cluster temperature and shape (via the 6 c dp^/ Ry factor). 

The variation of / gas with f^o, for an open Universe (A = 0) and a flat Universe (fio + A = 1) is plotted on Fig. 5 for 
the X-ray best fit parameters of MS1054-03 (kT = 12.3 keV, p = 0.96, 6 C = 0.89'). 
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Table 1 

Isothermal /3-model fit results 

#0 r c i r c2 (3 x Q y B 

(ct/s/arcmin 2 ) (arcmin) (arcmin) (J2000) (J2000) (ct/s/arcmin 2 ) 

Td (5.4 ±0.6) IO" 3 0.89l^| 0.96^ 22 (3.07 ±0.02) IO" 3 

2-d (6.0 ± 0.4) 10~ 3 0.85 ±0.22 0.60 ±0.13 0.73 ±0.18 10 /l 56 m 58.6 a -03 d 37 m 36 s (3.00 ± 0.03) IO" 3 



Table 2 
Mass estimates 



Method 


Radius 


M gas 


M tot 




(Mpc) 


(10 14 h" 5/2 M Q ) 


(10 15 M ) 


1- D BM 

2- D BM 


1. 
1. 


1.11 ±0.05 
1.27 


1 - 1 -0.2 n 50 

0.89 


1- D BM 

2- D BM 


1.65 
1.65 


1.9 ±0.3 
2.5 


^•U_ .4 n 50 

1.6 


VT 


1.65 




1.6 h^ 1 
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Fig. 2. — Restored X-ray image from the wavelet analysis of the HRI raw image. The iso-intensity contours of the reconstructed main 
cluster and of the compact sub-substructure (in the west) are overlayed on the optical image (D98). The isocontours are linearly spaced by 
6.1 10~ 4 ct/s/arcmin 2 , the value of the lowest contour. 
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Fig. 3. — Surface brightness profile of MS1054. The crosses show the data with la error bars. The full line is the best fit 1-D /3-model 
The dotted curves are the profiles of the best fit 2-D elliptical model ±2<r errors on the central surface brightness (j3 = 0.73). 
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Radius (Mpc) 

Fig. 4. — Total and gas mass profiles of MS1054. The error bars correspond to the error on the imaging data. Full lines and dashed lines 
are the results of the 1-D and 2-D isothermal /3-model fit respectively. Dashcd-dotted lines: virial mass from numerical simulations at this 
redshift (Evrard et al. 1996). The heavy lines are for kT = 12.3 keV. The thin lines correspond to the lower (kT = 10.1 keV) and upper limit 
(kT = 15.4 keV) of the ASCA temperature measurement. 
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Fig. 5. — Variation of the estimated gas mass fraction / gaa of MS1054-03 with the assumed cosmology. / gaa is normalized to the value 
obtained for Qq = 1. Full line: variation with Qq for an open Universe (A = 0). Dotted line: flat A dominated Universe. The X-ray 
parameters are fixed to their best fit values (kT = 12.3 keV,/3 = 0.96, 6 C = 0.89'). 
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